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' The phase structure of the 3D 5C/(2)-Higgs model, the dimensionally reduced 

effective theory for the electroweak model at finite temperature, is analysed 
■ on the lattice using a variant of the linear ^-expansion. We develop a sys- 

tematic variational cumulant expansion for general application to the study 
| of gauge invariant operators in 3D gauge-Higgs models, with emphasis on 

0^ ' the symmetric phase. In particular, the technique is not restricted to finite 

lattice volumes, and application to the fundamental 3D 5[/(2)-Higgs model 
^p -1 ' allows the discontinuity of certain observables across the first-order transition 

Oh. to be observed directly for small 4D Higgs masses. The resulting phase struc- 

t—I ! ture agrees well with Monte Carlo simulations for small Higgs masses, but, at 

least to the order calculated, the technique is less sensitive to the expected 
evolution of the transition to a crossover for Higgs masses above 80 GeV. 
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I. INTRODUCTION 



Recently, considerable effort has been invested in the study of the phase structure of the 
electroweak model and various extensions at finite temperature. In particular, the properties 
of the first-order phase transition are of great cosmological significance, in order to determine 
whether the generation of baryon asymmetry is viable at such a transition. 

The study of gauge theories, such as the electroweak model, at finite temperature is 
plagued by the fact that combined with the need for non-perturbative techniques to study 
the confining symmetric phase, infra-red effects also lead to a large effective expansion 
parameter in the broken phase where perturbation theory appears naively applicable. This 
is particularly true for moderately large Higgs masses, m# ~ mw- A considerable advance 
for the study of static, thermodynamic, properties has been the systematic development of 
dimensional reduction techniques [lj] allowing the mapping, via matching of Green functions, 
of the full 4D theory onto a super-renormalizable 3D theory corresponding to the Matsubara 
zero-modes. For the electroweak theory in the transition region this procedure may be 
carried out perturbatively as the gauge coupling is small, while infra-red problems are also 
absent as such effects are encoded in the dynamics of the effective 3D theory. If one ignores 
the unimportant U(l) factor ||, the effective theory is a super-renormalizable S77(2)-Higgs 
model. This super- renormalizability, and also the non-triviality of the Higgs sector in 3D, 
allow a more straightforward approach to the continuum limit than in conventional 4D 
studies |3|-|S| . Removal of thermally massive fermions, and smaller 3D lattices have allowed 
detailed Monte Carlo studies of the 3D theory, the conclusion being that the first-order 
transition which grows increasingly weaker as the Higgs mass increases, actually ceases to be 
first order and is presumably an analytic crossover for m# > 80 GeV ||9|-11]. This is plausible 



as in this system there is no gauge invariant order parameter which can physically distinguish 



the two phases [12]. Indeed recent analysis of the spectrum for large Higgs masses |]T3"|Jbl 
has indicated that the structure is very similar on both sides of the transition/crossover 
region. One can intuitively map three massive vector bosons and one scalar of the Higgs 
phase smoothly onto the three vector and one scalar bound state resonances of the confining 
phase. However, actual identification of the spectrum requires more care due to operator 
mixing and the presence of low lying excited states. 

Notwithstanding the success of Monte Carlo simulations, it is also clear that many of 
the advantages of dimensional reduction for the study of finite temperature gauge theories 
are more generally applicable. Indeed, renormalization group [ 15,16] , and Schwinger- Dyson 
|I7J , techniques were among the first to suggest loss of the first -order transition, stimulating 
much of the recent lattice activity. Furthermore, analysis of the spectrum, and properties 
of the crossover very close to the endpoint of the first-order transition line, may well prove 
difficult to study via Monte Carlo techniques due to the requirement for large volume lattices 
in order to overcome finite size effects. This is especially true if the transition line ends at 
a second-order endpoint. 

In this paper we shall investigate the fundamental 3D 5't/(2)-Higgs model on the lattice 
using an analytic technique which is a variant of the linear ^-expansion (LDE). The ap- 
proach is developed as a systematic method applicable to the calculation of gauge invariant 
expectation values in gauge-Higgs theories on lattices of arbitrary (even infinite) volume. 
We concentrate in this instance on the fundamental Si7(2)-Higgs model as the dimensionally 
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reduced effective theory for the electroweak model and consider the calculation of average 
plaquette and hopping term expectation values in order to extract the phase structure. Some 
preliminary results of this work were presented in |[18||. A similar methodology has been used 



previously for studying the average plaquette energy in 4D lattice gauge theory [|T^-|25|, the 
phase structure in a mixed fundamental/adjoint 57/(2) model pf] , Higgs models p7|-p9|, 
and calculation of the scalar glueball mass in pure SU(2) gauge theory [j3"0fl . The results 
have generally agreed well with Monte Carlo data where available. 

The layout of the paper is as follows. In section 2 we briefly review the relevant relations 
between the dimensionally reduced 3D theory and the electroweak 4D parameters, which 
allow us to present the results in terms of 4D Higgs mass and temperature variables. In 
Section 3 we review the LDE approach, and describe the application to gauge-Higgs systems, 
and the fundamental 5£/(2)-Higgs model in particular. The calculational techniques are 
discussed in section 4, with more technical details relegated to two appendices. Results 
for relevant observables in the transition region are presented, with higher order cumulants 
being used to help locate the critical line precisely; the results agree remarkably well with 
Monte Carlo estimates. Some concluding remarks are presented in Section 5. 



II. DIMENSIONAL REDUCTION OF THE ELECTROWEAK MODEL AT 

FINITE TEMPERATURE 

A full analysis of dimensional reduction of various 4D theories including the standard 
model and MSSM, involving detailed matching of Green functions has been presented in 
ffl. Here we shall simply summarise the procedure and present the relations between the 
parameters for later use. 

The reduction procedure consists of two stages. In the first, high momentum modes 
of the gauge field, and the thermally massive fermions are integrated out, while retaining 
the temporal component of the gauge field A as an additional adjoint Higgs. The mass 
of this field is given by the Debye screening mass as O(gT), and this field may also be 
integrated out as a second step. One then matches the resulting Green functions with those 
for a super-renormalizable 3D theory which, neglecting the U(l) factor, is an £7/(2)-Higgs 
theory, 

C = \g%G% + (A$) f (A$) + mj&Q + A 3 ($ t <l>) 2 , (1) 

where is the 57/(2) field strength, and <3> an 577(2) doublet. More recently, study of 
the full 57/(2) x U(l) theory has indicated @ that the phase structure is not qualitatively 
changed, thus for convenience we consider only the case where the Weinberg angle is set to 
zero. It is an important feature of the Green function matching procedure that it does not 
introduce additional non-local terms which generically arise in a low energy effective theory 
obtained by explicitly integrating out massive modes, and may give dominant contributions. 
It has been shown Jl| that the Green functions of ([]]) approximate those for the full theory in 
the infra-red up to corrections of the form 5G/G ~ 0(g 3 ) where g is the 4D gauge coupling. 
This is an acceptable approximation as the phase transition region is expected to have a 
weakly coupled gauge sector. 
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The 3D gauge coupling g%, Higgs mass m\ and self coupling A3 will depend on the tem- 
perature and the underlying 4D parameters. In this theory all parameters are dimensionful 
and one may fix the scale with the gauge coupling g\. The phase structure and 4D temper- 
ature dependence are then determined by the two dimensionless ratios [[J x = A 3 /g 2 and 
y = m\(g1) / ' g\. However, the 3D theory is an effective theory for a large class of 4D theories 
with differing scalar and fermionic field content, and the mapping needs to be established for 
each individually. As discussed above, this may be performed perturbatively as calculations 
only involve massive modes, and are thus free of infra-red problems. For the standard model 
the relevant mapping has been determined in JT]. For simplicity, we consider the mapping 
to the effective 4D parameters m* H and T* of the 4D S'?7(2)-Higgs model |]3T|,|9[| which are 
convenient for representing the results, and differ by only a few percent from the pole mass 
and temperature. Using tree-level relations [pif , the temperature is given purely by the 
gauge coupling, g\ = 0.44015T*, while the dimensionless ratios are 



x = -0.00550 + 0.12622/i 2 

(M* ) 2 

y = 0.39818 + 0.15545/i 2 - 0.00190/i 4 - 2.58088 K ,J I > . 

(T*y 



(2) 
(3) 



where h = m^/80.6 GeV is a dimensionless Higgs mass. As there is this direct connection 
between (m* H , T*) and gf , x, and y we can convert analysis of the phase structure in the x, y 
plane to the physically more meaningful (m* H ,T*) variables. 

Focusing now on the 3D theory, the fact that the infrared details are still encoded 
in the dynamics of this theory means that a non-perturbative approach is required. The 
super-renormalizability of (P now ensures that one can calculate perturbatively exact RG 
trajectories. Only the mass parameter receives contributions at two-loop order, and there are 
no three- or higher-loop contributions. We adopt a lattice regularization, and the standard 
Wilson lattice analogue of the action (|l|) with lattice constant a may be written in the form 



S 



p z 



p + 2& £ Pip,™* -£[/>•+ MpI - if 



(4) 



where we choose to represent the Higgs doublet in the form $ = pjVj, with p GR + and V an 
element of the fundamental representation of SU(2). We then make a gauge transformation 
Uij — > ViUijVj to rotate the phase of the Higgs field in the interaction term to the identity 
at each site. Thus the interaction term becomes pipjTiUij. While this choice may not 
appear desirable in the symmetric phase, it is advantageous for calculational purposes, and 
in previous lattice studies P7| has led to good results, so we shall also adopt it here. In 
the continuum limit, which may be determined by the RG trajectory noted above, the 
dimensionless parameters /3, flh, j3 r are related to the parameters of the continuum theory as 



follows [31 



x 



+ 



16tt 2 



-3 



/51 

— + 9x 
V16 




(5) 



5.2x 



(6) 
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where £ = 3.17591 and ( = 0.09 follow from perturbative analyses. For different values of 
the lattice constant a these relations define a constant physics curve, or RG trajectory, for 
fixed g^,x,y in the space j3, j3h, j3 r . 

With the preceding results, we may now analyse the phase structure of the lattice SU (2)- 
Higgs model in 3D and subsequently convert the results back to physical parameters. The 
details of this procedure will be discussed in Section 4, while in the next section we review 
details of the LDE, and discuss its application to the SU (2)-Higgs model. 



III. LDE ANALYSIS OF THE S77(2)-HIGGS MODEL 



The linear 5-expansion is an analytic approach to calculations in field theory which makes 
use of an artificial expansion not dependent on the existence of small coupling constant 
for its validity. Nonetheless, calculations are generally no more complex than standard 
perturbative Feynman diagrams. An essential feature of this approach is the ability to 
optimize convergence of the series by fixing additional parameters appearing in the extended 
action. 

The LDE action has the form 



S 5 = (1-5)S (J)+SS, 



(7) 



where So contains some dependence on a variational parameter J, and S$ — > S the theory 
under consideration, independent of J, when 5 = 1. The generating functional for Green 
functions may then be expanded to an appropriate order in 6, which is then set to unity. 

As the power series is only calculated to a finite order, it retains some dependence on J 
which would be absent in a full summation. A well motivated criterion for fixing J is the 
principle of minimal sensitivity (PMS) ||32|| , whereby J is chosen at a local extremum of a 
physical quantity, i.e. if X N (J) denotes the iV'th approximation to X, then we impose 



dX N (J) 
dJ 



0. 



(8) 



This or a similar criterion is intrinsic to the success of the LDE, providing the nonper- 
turbative dependence on the coupling constant, and has been shown to induce convergence 
of the series in 0- and 1-D field theories |33,33] where the perturbative series are asymptotic 
and eventually diverge factorially. 

The technique has also been applied with success to lattice gauge theories, where different 
choices of Sq are appropriate in the weak and strong coupling regimes. In the strong coupling 
regime an approach using a maximal tree of plaquettes [Pljpll has proved successful in 
describing the strong coupling behaviour of the average plaquette energy Ep. The weak 
coupling regime has also been considered using a quadratic So |]2"2"| , |2"B"|] and using a single link 
These techniques have since been applied at finite temperature | 20fl , to the 577(2) 
mixed S77(2) 



So in 



mass gap pO 
in 



50(3) phase structure [P6[ , and also to lattice Higgs models 



P7H2 



In the present context, the variational cumulant expansion approach using a single link 
So seems most appropriate. At this stage we retain the full gauge Higgs coupling and choose 
a background action of the form (cf. (§)) 



5 



So = E( J + (3 h PiPi) l ~ftUi - E [a 2 + &(p 2 - 1) 2 ] , (9) 

I i 

thereby introducing a variational parameter J. It should be noted that this background is 
not explicitly gauge invariant due to the single link action. In the present context this is not 
of concern since we shall only calculate the expectation value of gauge invariant operators. 
However, if one were to calculate gauge non-invariant expectations the correct approach 
20| is to note that a priori J is link dependent and may have an arbitrary SU(2) phase 



which should be summed over. Once this procedure is carried out one finds that gauge 
non-invariant expectation values vanish as expected^. Using the methodology just described 
we evaluate expectation values as follows: 

(X) = ^ f [dU}{dp\Xe 5{s - So) e s \ (10) 
where the partition function is given by 

Z = e~ w = J[dU][dp\e 5 ^- s ^e So , (11) 

in which [dU] denotes the standard Haar measure over SU(2), and [dp] = TiidPiPf- An N th 
order approximant to (X) in the LDE then takes the form, 

N 5 n 

( X )n =T,ZJ(X(S-S°) n )c> (12) 

n=0 n - 

where (X) c denotes the connected expectation value, or cumulant, of X in the So back- 
ground, e.g. (AB) C = (AB) Q — (A) (B) . It is the expansion of the partition function in the 
denominator of the expectation value to the appropriate order in 5 that naturally subtracts 
the disconnected components. This is crucial, as one now observes that the physical expec- 
tation value reduces to a sum of connected expectation values for which one can introduce a 
convenient diagrammatic notation in order to keep track of the independent contributions. 
Furthermore, connectedness implies that one can effectively take the lattice volume to infin- 
ity as at any finite order only a finite number of distinct diagrams will contribute. Summing 
all such diagrams over an infinite lattice will then simply lead to coefficients proportional to 
the number of sites N in the lattice, and thus expectation values normalised as (X) /N will 
be finite and calculable, independent of the finiteness or otherwise of N. 

Therefore, for practical purposes the calculation reduces to the enumeration and calcu- 
lation of expectation values of connected diagrams in the S background, i.e. we require^] 



It is interesting to note that this approach does in principle allow the possibility of considering 
gauge fixed operators, and this may be of interest in comparing gauge invariant lattice results for 
the spectrum with gauge fixed perturbative calculations [f35f-|37]]. 

2 Note that if convenient one may reduce the number of diagrams required by calculating (X(S n )) 
in an Sq(1 — 5) background, and then subsequently expand the result to the appropriate order in 
5. 
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<X> =1- ({dU}{dp]Xe s \ (13) 
z J 

with 

Z = f[dU][dp}e s °, (14) 



Actual calculations are simplified by recognising that the terms in the S*?7(2)-Higgs 
action involve only real constants, the p, n or SU(2) characters since TrU = Xi/2{U) for U an 
element of the fundamental representation of SU (2). The free action Sq also consists only of 
constant factors and fundamental SU (2) characters. For example, the Sq partition function 
takes the form, 



Zo = /[II ^Pf\ J til exp fe( J + /Wi)Xi/2(£%) - E VhOSt, Pi)) , (15) 



where Vh = p 2 + (3 r {p 2 — \) 2 is the Higgs potential. The expectation value integrals are 
most easily evaluated making use of a character expansion for the exponentiated Sq. For 



the fundamental representation of SU (2) we have (see e.g. |38 



eiJxi/adO = ^ £ (r + l)/ r+1 (J) Xr/2 ([/), (16) 

^ r=0 

where we have made use of the relation I p (x) — I p+ 2(x) = 2{p + l)I p+ i(x)/x for modified 
Bessel functions. Then via use of the orthogonality relation for characters the group inte- 
gration becomes trivial. In certain higher order cases expectation values do not reduce to 
simple products of link characters and a rather more general technique is then required (see 
Appendix A for more details). In order to make the calculations tractable we will make the 
assumption that the Higgs modulus varies relatively slowly over the lattice. More precisely, 
we assume that piPj = p 2 + 5 2 piApi, where Api = pj — pi with i and j nearest neighbour 
sites. This is motivated by the success of the approach, ignoring the correction altogether, in 
spin-glass systems |39|, which are very similar to the (3 — > limit of the S"?7(2)-Higgs lattice 
theory. The validity of this assumption for large /3 is less clear. However, given our ansatz, 
at higher order in 5 the correction terms enter in the calculation in a self-consistent way. 
Consequently, up to 0(S 2 ), once the group integration has been performed via character ex- 
pansion techniques, the remaining integrals over the modulus of the Higgs field generically 
take the form (for each link) 




l/d 

(17) 



which may be calculated numerically (see Appendix B for an analytical study, and an ex- 
pansion in (3h). All the required expectation values may then be expressed in terms of the 
ratios B n r = A?/A° . 
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IV. EXPECTATION VALUE CALCULATIONS AND THE PHASE STRUCTURE 



In considering expectation values it is convenient to introduce a graphical notation for 
the standard operators appearing in the action. We define: 

° = E^™p (18) 

— = Y.\piPi^ u a ( 19 ) 
i z 

-=E^i- (2°) 

The gauge invariant quasi-order parameters, the average plaquette, and the hopping term, 
may then be conveniently represented in the form 

(Ep) = ^ (□) (21) 

(El) = ± (— ) , (22) 

where the expectation value in this case corresponds to the full action and, working with an 
TV-site (i-dimensional lattice, the number of links and plaquettes is given by N t = Nd, and 
N p = d(d — l)N/2, respectively. 

Disconnected contributions to these observables in the So background are calculated in 
the manner introduced in Section 3, and elaborated further in Appendix A. The required 
expectation values, along with the associated multiplicities, of the relevant diagrams are also 
tabulated in Appendix A for reference. In order to present reasonably concise expressions 
we refer to the relevant connected contribution from diagram di as Cj = TOj (Dj) c where mi 
is the multipicity and (Di) c is the connected expectation value obtained by subtracting the 
disconnected components from (Di) in the standard manner. 

We consider first the pure gauge sector of the theory. It is clear that in the limit (3h — > 
the partition function reduces to that of pure SU (2) gauge theory multiplied by an overall 
factor determined by the integration over the Higgs modulus, which factors out of expecta- 
tion values. We now consider the average plaquette for non-zero (3h in the symmetric phase. 
To first order in 5 the expectation value is given by 

(E P ) = i-«D) + 5 (n(S - S )) c + 0(5")). (23) 

Recalling that only connected diagrams contribute, the contribution to 0(5) is given by 

(Ep), = d + P(C 6 + C 7 ) - JC 8 . (24) 

The variational parameter must now be fixed by using the PMS criterion described 
earlier. Considering for a moment a general observable (0(/3, /3h, Pr, J)), the procedure used 
throughout is to examine (O) at fixed (3,(3h,(3 r as a function of J, fixing the parameter at 
the "flattest" local extremum J = J in accordance with PMS. As a result, the expectation 
value has the functional dependence (O) = (0({3,(3h,(3 r , J(P, Pt» Pr))/ an d, due to 



8 



FIGURES 




FIG. 1. The plaquette energy E p as a function of for (f3h = 0, f3 r = 0) and also {j3h = 0.2, 
Pr = 0.5), a smooth transition between PMS branches occurring at (5 ~ 2.44 and 2.37 respectively. 

discontinuous jumps between the PMS extrema J, may exhibit non-analyticities even when 
(0(/3, fit, Pr, J)) is an analytic function of its parameters. Thus one associates the non- 
analyticity of discontinuous observables across a first order transition with jumps in the 
PMS solution J. Therefore the appearance or disappearance of a particular J as one moves 
over the parameter space is associated with a phase transition of the system if the resulting 
observable exhibits a discontinuity. This interpretation is justified by the experience with 
this technique in lattice gauge theory and in particular the study of the SU{2) — SO (3) 
phase diagram p6 |, where the results agreed very well with the expectations from Monte 
Carlo simulations. 

Returning now to the average plaquette expectation, we plot the results in Fig. 1 as a 
function of (3 for f3h = (pure gauge theory) and (3h = 0.2 (well into the symmetric phase). 
One observes that the results are remarkably similar, indicating that the presence of the 
Higgs sector has little effect. This conclusion is consistent with explicit analysis of the (Ep) 
expectation values (see Appendix B), which suggests that the first Higgs sector corrections 
in the symmetric phase are generically of 0((/?/ l /2) 4 ). 

The structure of the free action So indicates that the limit J = corresponds to the 
standard strong coupling expansion, and indeed a PMS solution exists for all (3 at J = 0. 
A second PMS solution first appears at (3 ~ 2.4, then flattens out to become the clear 
PMS point as f3 increases. This transition to a "weak coupling" branch occurs, however, 
at almost precisely the same value of (Ep), and with no distinctive change in behaviour, 
implying that the transition is actually a crossover. This deconfinement crossover is indeed 
consistent with expectations from lattice Monte Carlo studies in AD 0H|J, as is the slight 
shift of the crossover to lower values of (3 as fih increases. 

We now turn now to consideration of the hopping term which, for small Higgs masses, 
serves as a useful gauge invariant operator which is discontinuous across the first-order 
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transition. In order to analyse the behaviour as precisely as possible we take the calculation 
toO(5 3 ), 



1 / 3 5 n \ 
(El) = w £ - (—(^ - S T) C + 0(5') . (25) 

^1 \n=0 n - J 

The relevant connected diagrams and their expectation values are enumerated in Appendix 
A. The result is given by 

(E L ) =C 1 + 6{2(d-l)flC i -JC 5 ) 

+y (V^(Ci + C n + C 12 ) - (3J(d - 1)(C 13 + C 14 ) + J 2 C^ 

A3 / A _ 1 36 r/ - 1 43 

\ Z i=27 Z i=37 

+3/?J 2 ^(C 44 + C 45 + C 46 ) - J 3 C 47 ^) + 0(5 4 ). (26) 

Considering first the expectation value at 0(5) only, we indeed find a transition between 
PMS branches in the region of (3h m 0.31 — 0.38 consistent with Monte Carlo expectations 
for the first-order transition. However, there is a region in parameter space in the transition 
region between these two branches for which no PMS extrema exist. This is not unexpected, 
as the expectation value undergoes a discontinous shift. However, it does mean that a precise 
determination of the critical point will prove difficult working directly from (El). Analysing 
a range of f3 values, the results are presented in Fig. 2 for a 4D Higgs mass of 70 GeV, where 
this constraint determines (3 r via (H) and @. Since the approximation used in evaluating 
the integral over the Higgs modulus is expected to be most accurate in the j3 — > limit, we 
restrict ourselves to relatively small values of (3 in this particular calculation. 
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FIG. 2. We plot < El >i versus (5^ for a 4D Higgs mass m* H = 70 GeV, and three values of f3. 
The vertical line is simply a rough guide to the position of the transition region (3^ = 0.31 — 0.35. 



10 



The results presented in Fig. 2 are extended to 0(5 3 ) in Fig. 3. One observes that the 
transition region is now more clearly defined, with the PMS branches persisting closer to the 
transition point itself. At this order one has additional terms from diagrams at 0(5°) (dl), 
and 0(5 1 ) (d4 and d5), correcting the slowly varying Higgs approximation in the evaluation 
of A™. However, we note that in the transition region the behaviour of (El) is dominated by 
long wavelength modes, and thus the approximation used should be appropriate. A sample 
analytic calculation of the correction term is presented in Appendix B, and provided one 
assumes that the modulus of the Higgs field varies by no more than 10% over a lattice 
spacing, one finds that the correction terms are negligibly small when compared with the 
original contribution at the appropriate order in 5. 






0.2 0.3 0.4 0.5 

FIG. 3. We plot {El) 3 versus (3h for a AD Higgs mass m* H = 70 GeV, and three values of j3. 
The vertical line is at the same value of (5h as in Fig. 2. 

In order to determine the transition point more precisely, we now consider the second- 
order cumulant of El, defined by C(El) = (El) — (El) 2 . To 0(5) this expectation may be 
represented as, 

= C m + S \P^-{C m + C50 + C 51 ) - J(C m + C 53 )j + 0(S 2 ). (28) 

The cumulant exhibits a PMS solution for small /3h, but as /3h is increased this solution is 
lost at /3hc, above which no PMS extrema exist. Accordingly, we associate this point with 
the phase transition point. Tracking this line for increasing j3 we obtain an estimate for the 
transition line along the RG trajectories in the phase diagram for 4D Higgs masses of 60 
and 70 GeV. The results are shown in Fig. 4 and compared with known Monte Carlo points 
in the diagram from the simulations presented in []3l], | |4"0]| , and Lee- Yang zeros analysis in 



ill . Our procedure, even at this low order, produces close agreement with the Monte Carlo 



results. 
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FIG. 4. We plot the critical curves in the (f3h,P) phase diagram for Higgs masses of 60 and 70 
GeV. The critical values of (5h were determined to be those points at which a PMS extremum for 
the cumulant C(El)i was just lost; the values corresponding to a discontinuity of the cumulant. 
The results are compared with the Monte Carlo results of [31] (squares), and Q] and JT1J (crosses) 
for m* H = 60 GeV. 

Extrapolating to the continuum limit along the RG trajectory, we obtain critical cou- 
plings for each Higgs mass. Using (^) and (|6[) we can also determine the critical temperature 
and perform a similar extrapolation to the continuum. However, as one observes from Fig. 4, 
the continuum limit of /3h c is approximately 3% above the Monte Carlo estimates. Although 
this is entirely satisfactory for the order of the calculation, the sensitivity of the critical 
temperature to the value of (3h c results in estimates for the critical temperature being sig- 
nificantly lower than one would expect. In the low f3 region, however, where the results 
at this order are expected to be most accurate, the critical temperature agrees with Monte 
Carlo results |3IJ up to 20-25%. Note also that while the first-order cumulant appears to 
over-estimate the critical coupling j3h c , the transition region determined by analysis of (El) 
gives an under-estimate, for the relevant (3 range, of 6-8%. This discrepency provides a 
reasonable estimate of the systematic errors of these results. 

Increasing the Higgs mass beyond the regime where the first-order transition is expected 
to end, we find that the results change qualitatively only in the sharpness of the transition. 
In Fig. 5 the observable (El) 3 is plotted at (3 = 3 for a AD Higgs mass of 120 GeV. One 
observes that the PMS solutions do not persist as far into the transition region in comparison 
with the results at m* H = 70 GeV, and that the discontinuity itself is slightly smaller in 
magnitude. Although it is not possible from this result alone to characterise the transition 
as an analytic crossover rather than a first-order transition, it is certainly consistent with 
the former interpretation, which is suggested by previous studies. 
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FIG. 5. We plot < E L > 3 versus fa for a 4D Higgs mass m* H = 120 GeV at (5 = 3. The vertical 
line indicates the transition/crossover region. 

The results are also consistent with the fact that Monte Carlo simulations indicate that 
the crossover is still extremely sharp for a Higgs mass in this region. Fig. 6, in which the 
calculation is repeated at m* H = 140 GeV, shows a further reduction in the sharpness of 
the transition, as one observes from the increase in the parameter range for which no PMS 
solution exists. 
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FIG. 6. We plot < El >3 versus fa for a 4D Higgs mass m* H = 140 GeV at (3 = 3. The vertical 
line indicates the transition/crossover region. 

Increasing the Higgs mass beyond ~ 150 GeV leads to a loss of the PMS solution in the 
broken phase. This is most likely to be a finite order effect with little direct relevance for the 
true phase structure, and for this reason we shall ignore such large Higgs masses. In studying 
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the second-order cumulant, we find that the qualitative behaviour is again unchanged until 
the Higgs mass becomes very large. In the region of interest, 80 — 120 GeV, the loss of a 
PMS solution at a critical value of f3h may presumably be associated with the sharpness of 
the crossover in this regime. 



V. DISCUSSION 

In this paper we have made use of a variant of the linear (^-expansion to examine the 
phase structure of the 5'C/(2)-Higgs model in 3D. Using the previously derived matching 
conditions between the parameters of this model and a suitable approximation of the full 4D 
finite temperature electroweak theory we were able to extract the phase structure, consistent 
with lattice Monte Carlo simulations for small Higgs masses. Increasing the Higgs mass 
beyond the point where the first-order transition is expected to be replaced by a smooth 
crossover, the discontinuities in the observable (El) persist, but become less pronounced, as 
m* H increases. A higher order calculation in S may be necessary in order to observe a more 
significant qualitative change in the behaviour of the relevant observables. 

Nevertheless, this work does suggest a number of immediate extensions which are cur- 
rently under consideration. The formalism developed here, in combination with the approach 



used in P0| , allows an analytic study of the spectrum in the symmetric phase, in particular 
near the critical endpoint, where the lack of any finite volume restriction makes the regime 
dominated by large correlation lengths quite accessible without detailed finite volume scaling 
analysis. Recent work on the spectrum of the 3D theory for large Higgs masses suggests 
that in this respect the Higgs and confining phase have very similar properties [II . 



The technique also suggests the possibility of studying the adjoint 5*^7 (2)-Higgs model, as 
the dimensionally reduced effective theory for QCD f4T|-f43|| . A study of the phase structure 
of the mixed fundamental and adjoint SU(2) gauge theory f26|| , which is the lowest order 



correction to pure SU(2) gauge theory from a weak coupling expansion in fih of the adjoint 
Higgs model ^TJ, has already been successfully carried out using similar techniques. 

However, given the remarkable level of precision that Monte Carlo techniques have been 
able to attain in this area, it is perhaps more appropriate to pursue the study of problems 
which are rather less accessible to such simulations. For instance, the possibility exists for 
studying the full AD finite temperature theory directly J2U|. Furthermore, study of the 
[/(l)-Higgs model, with the inclusion of a chemical potential in the context of superconduc- 
tivity, appears feasible whereas, in contrast, the additional phase introduced by the chemical 
potential renders Monte Carlo simulations for this system very problematic. 
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APPENDIX A: CALCULATIONAL TECHNIQUES AND EXPECTATION 

VALUES 



In this appendix we briefly discuss certain aspects of the techniques used for evaluation 
of the expectation values within the LDE. All the required expectations may be reduced to 
the form 

<X> =i- f[dU][dp]Xe s °, (Al) 

where in all cases X corresponds to a product of real constants, factors of p, and 577(2) 
characters Xr{U) = (TrU) r . In addition 5*0, up to terms in p, also consists of fundamental 
SU (2) characters. Thus the exponential may be represented as a character expansion 

e* J *' aW ^£W-0Xr/a(E0, (A2) 

where I r (J) are modified Bessel fuctions. 

Almost all the required expectation values X may be converted into simple character 
factors for each link using a combination of Clebsch-Gordan expansions, and the following 
relations, which follow easily using the invariance of the Haar measure and Sq under a 
unitary transformation ||19||: 



< Xr/2{U 2n ) > = (r + l) 1 - 2 " < Xr /2(U) > 2 " . (A3) 

Note that these relations also hold when there are additional single-link characters of U 
present in the expectation. With these techniques most of the required expectations are 
reduced to a product of constants and simple character factors. The group integration then 
becomes trivial via use of the orthogonality relation 

[du\xr{uv)xAu- x w) = -£—gr' Xr {yw). (A4) 

The remaining integration over p may then be evaluated numerically using ([H|). As an 
example, consider evaluation of the following expectation value: 

Wo = ~ (p 2 1 Xi/2(U 1 U 2 uM)Xi/2(U 1 )) o (A5) 
= i (p\xi/2{U 2 )Xi,2{Ui)xi/2{Ul){l + Xi(U 1 ))) Q (A6) 
= ^(5?) 3 (5o 1 + 5 2 1 ). (A7) 



In the second line we have made use of (A"3), and a Clebsch-Gordan expansion for (xi/2(^7)) 2 - 
However, there are certain diagrams appearing at third order in 5 which cannot be 
reduced to simple character factors for each link due to the degree of connectedness of 
the plaquettes. For example, see diagrams d20, d35 and d36. In order to calculate the 
required expectations in this case we use the projection operator technique elaborated in 
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which allows the direct evaluation of un-traced products of a single 
for more details. 



30 



the Appendix of |30 
link. We refer the reader to 

For each diagram required for the analysis in Section 4 the disconnected expectation 
values ((-Dj)) calculated in the manner described above, and also the corresponding mul- 
tiplicities (mj), are now presented in the following table. For convenience we define the 
constants d\ = 2d — 3, d,2 = 2d — 4, = d — 1, and d^ = 2d — 1. 
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d# £i rrii (A)q 

d21 [2] | 4 l^ + M^nBV + Bl)) 

d22 I I 1 1 24rfi ^(S?) 5 (l + S§) 2 

d23 []□ 4dx £(a?) 8 (25° + 5°) 

d24 Q| 4 £(iJ?)S(25° + a°) 

d25 |Q 12 £(a?) 2 (l + B°) 2 

d26 HI 1 1(25° + 5°) 

d27 [3]J 4 ^(^) 3 (S 1 + 5 2 1 ) + ^(S 3 ) 3 (S 2 1 + 5 4 1 ) 

d28 [2l~| I 36rf i 2 ^(fi 1 ) 2 (fi 1 + fi 2 1 )(fi 1 + ^(5 2 ) 3 (fi? + fi 3 )) 

d29 Q|]J 36d! ^(S^^Si + ^^Sy + BSJCSi + Si)) 

d30 rl^l 12rf i 2 V(5 1 °) 3 (S 1 + J B 2 1 + ^(S 2 °) 3 (S 1 + 25 2 1 + J Bi)) 
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d38 @J 12 ^(5o 1 5 1 + Jr(5 2 ) 2 (5 1 + 53 )(5 1 1 + fi 3 1 )) 

d39 m il 24rf i 2 4o(5?) 5 (l + J B 2 °)(2fi 1 1 + ^) 
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APPENDIX B: EXPANSION OF A%, AND HIGHER ORDER CORRECTIONS 

After the group integrations have been performed, all expectation values up to 0(5 2 ) 
reduce to a series of ratios of the integral ([Lj 



iKY - fV^ ((r + D ^ff )'. (Bl) 

for various choices of n and r depending on the diagram being evaluated. While the integral 
may easily be performed numerically, it is useful to consider an analytic series solution, 
particularly in the strong coupling regime (3 — > 0, and in the regime of a minimal gauge- 
Higgs coupling, f} h < 1. 

As was discussed in Section 4, the LDE approach naturally includes a solution at J = 
which corresponds to the strong coupling expansion appropriate for ^ < 1. We now wish 
to express ([Bl]) as an expansion in (3^. Expanding the modified Bessel function in a power 
series, we may represent the integral in the form 



oo 

2n+3-V H 



{A n r ) d = / dpp 2n+i e 







.fc=0 

oo 



d 



(B2) 



dpp 2n+3 e- v « £ c m (B3) 

m=0 



where 



Xk ~ k\(k + r + l)\ [—T- ) ' (B4) 

and in (|B3|) we have expressed the d th power of the series in A^ as a new series with terms 
c m which satisfy the recursion relation, 

i m 

^2{pd-m + p)X p c m ^ p , (B5) 
P =i 



mA 



with initial condition cq = Aq. This recursion relation may be solved by noting that the 
structure implies that c m is given by a sum of products of the Afc whose exponents in each 
monomial are given by an appropriate partition of m. The structure of the coefficients of 
these monomials is then easily deduced, and one finds that the terms c m have the general 
form 

Em 

c m = E A ^\ v=llP \\K P - (B6) 

{i p } 1 l p =i i P - [ d ~ 2^ P =i % v ) ■ p=i 

In the strong coupling limit with J = one may then readily perform the integration over 
p to obtain A™ as a power series in (3h- In the general case, after forming the product of the 
factors \ p p in (|BB|), we expand A^ in flBH| ) as a power series in @ h , giving c m as an explicit 
power series in /3 h . For compactness in the following results we define a t = SJ?Li (<?)%, where 
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the i q are determined by the relevant partition of m. With this definition, the result for c n 
becomes 



where 



{ip\ai=m} 



P(d, m, r) 



,jyd+2a 1 rd+2a 1 ( rd+2(Tl \ (fa 



(B7) 



s=0 



d!(r + i; 



(B8) 



( r \y-°°(d - <j )\ UT=i WKt + r + 1)0** ' 
Returning to the full integral, the p integration may now be performed using the relation 



/•OO i 

/ dpp 2x+3 e- VH ^ = Ux + l)\C{fa 2x), 
Jo 2 

where C may be written in terms of a parabolic cylinder function D n (x) as 

d-2/3,) 2 * f . i x (1 - OR \ 



Wr . 

which tends to unity as j3 r — > 0. Finally, the full integral is given by the power series 

. j x rd+2ai 



(B9) 



(BIO) 



/ / 

(^) d =E E 

™— n r„- i i \ z 



m=0 {ij,|(ji=m} 

rd+20"! 

X 2 S 

Z s=0 



Q(/3 r ,2n + 2s), 



(BH) 



where Q(/? r , 2a;) = (rr + l)\C(/3 r , 2x)/2. In the strong coupling limit (J = 0) this reduces to 

rd+2a 1 



(4?)s=E E ^(d,"»,r) 



Q(/3 r ,2n + rd+2(Ti; 



(B12) 



rra=0 {j p |(T 1 =m} V / 

The lowest order terms of this expansion are given by 



r ,nyf (P h Y d (n + rd/2 + l)\ 
[Ar)s -\2) 2(r\y 



1 + (t) , -^» + "'/ 2+2 ) +0 (t) 4 ' 



(B13) 



Thus, if we focus attention on the gauge sector and the lowest order contribution to the 
average plaquette energy, given by d2, we find that in the strong coupling region the f3h 
dependence is given by 



i + 



3 (Ph 
2 V 2 



O 



'fa 



(B14) 



Therefore, according to the expectation value (D 2 ), in this regime the leading dependence 
on j3h has the form (/^/2) 4 , which is a small correction throughout the symmetric phase 
where flh < 0.34, a result which is borne out by the numerical analysis presented in Fig. 1. 
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While the previous results are exact for expectation values up to 0(5 2 ), our restriction 
of the Higgs modulus in the interaction to the form p^pj = p? + 5 2 piApi implies that higher 
order calculations of observables will require additional correction terms for the lower order 
diagrams. As we have worked only to 0(5 3 ), it is only necessary to consider the first-order 
correction to the 0(5°) and 0(5) contributions to (El). 

Calculation of the correction is conveniently illustrated in the case of Zq. Expanding the 
exponential correction term we obtain, 



[dp] n / m 



i e 



|(J+/3fcP?)Xi/2(%) 



oo s.2n 

y — 

n=0 Z 



(B15) 



where we define [dp] = Hi d-PiPf ex P( — Vh)- On performing a character expansion of the 
exponential, Z may be expressed in the form 



[dpinx: E(r+i) 

lij n=0 rGZ+ 



I r+1 (J + p h pf) 5 



(PiPj - PIT / [d^;]Xr/2 



n/2 

E C jXj- 
=0,1/2 



(B16) 



The final sum begins with j = (j = 1/2) for even (odd) n, and the coefficients Cj are 
determined by the expansion of the product (xi/2(£^j)) n - This result is quite general and 
includes all higher order corrections. As mentioned earlier, we only require the first-order 
correction in the present case, corresponding to n = 1. To this order we have 



[dp] n 




1 + 5 2 {piPj - p\ 



'hiJ + PtP 2 



0{5 4 ). (B17) 



If we focus on the expression in square brackets, then since the ratio of Bessel functions has 
an upper bound of unity for large values of its argument, one observes that the term of 0(5 2 ) 
is indeed a small correction provided that p» ~ pj for i and j nearest neighbour sites. If 
one makes the assumption that pi and pj differ by at most 10%, which seems reasonable for 
expectation values dominated by long wavelength modes, then one finds that the correction 
is almost negligible despite appearing to be of 0(10%) in (P17| ). The reason is that, for 
non-zero f3 r , the integral ( |B17| ) is dominated by contributions from the region p < 10 — 15 
for which the ratio I2/I1 is small. Furthermore, the correction needs to be compared with 
the disconnected expectation value at the appropriate order in 5 rather than directly with 
its zeroth order contribution. Once one does this the correction is indeed found to be small 
having negligible effect on the results, as was discussed in Section 4. 

This result may be made more explicit by extracting terms up to 0(5 2 ) from the product, 
which take the form 1 + 5 2 J2<ij>(PiPj ~ Pi)h(30i)/Ii(xi), where x« = J + Phpf, 



S 2 



1 + T2 E z hj + 0(6* 

^0 <ij> 



The first factor is given by the standard zeroth-order contribution, 



Z = / [dp, 



(bis) 



(B19) 
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while the 0(5 2 ) correction also contains a sum of contributions of the form 

Z hj = f [d Pl dp,] { PlP] - p*) n (^) 4 jH- (B20) 

A further reduction in the relative magnitude of the correction would then follow from 
possible cancellations of terms in the sum over the links of the lattice, as the factor (piPj —pj) 
will vary in sign. 
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